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The 

RICIS 

Concept 


The University of Houston-Clear Lake established the Research Institute for 
Computing and Information systems in 1986 to encourage NASA Johnson Space 
Center and local industry to actively s upp ort research in the computing and 
information sciences. As part of this endeavor, UH-Clear Lake proposed a 
partnership with JSC to jointly define and manage an integrated program of research 
in advanced data processing technology needed for JSC’s main missions, including 
administrative, engineering and science responsibilities. JSC agreed and entered into 
a three-year cooperative agreement with UH-Clear Lake beginning in May, 1986, to 
jointly plan and execute such research through RICIS. Additionally, under 
Cooperative Agreement NCC 9-16, computing and educational facilities are shared 
by the two institutions to conduct the research. 

The mission of RICIS is to conduct, coordinate and disseminate research on 
computing and information systems among researchers, sponsors and users from 
UH-Clear Lake, N ASA/ JSC, and other research organizations. Within UH-Clear 
Lake, the mission is being implemented through interdisciplinary involvement of 
faculty and students from each of the four schools: Business, Education, Human 
Sciences and Humanities, and Natural and Applied Sciences. 

Other research organizationsare involved via the “gateway” concept, UH-Clear 
Lake establishes relat ionship s with other universities and research organizations, 
having common research interests, to provide additional sources of expertise to 
conduct needed research. 

A major role of RICIS is to find the best match of sponsors, researchers and 
research objectives to advance knowledge in the computing and inform|tipn: 
sciences. Working jointly with NASA/JSC, RICIS advises on research needs, 
recommends principals for conducting the research, provides technical and 
administrative support to coordinate the research, and integrates technical results 
into the cooperative goals of UH-Clear Lake and NASA/ JSC. 
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Abstract 

Periodic maneuvering control is developed for asymptotic moment um manage- 
ment of control moment gyros used as primary actuating devices for the Space 
Station. The proposed controller utilizes the concepts of quaternion feedback con- 
trol and periodic-disturbance accommodation to achieve oscillations about the con- 
stant torque equilibrium attitude, while minimizing the control effort required. 
Three-axis coupled equations of motion, written in terms of quaternions, are de- 
rived for roll/yaw controller design and stability analysis. It is shown that the 
quaternion feedback controller is very robust for a wide range of pitch angles. It is 
also shown that the proposed controller tunes the open-loop unstable vehicle to a 
stable oscillatory motion which minimizes the control effort needed for steady-state 
operations. 
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Introduction 


The Space Station will employ CMGs (control moment gyros) as primary actu- 
ating devices during normal flight mode operation. Gravity-gradient torques will 
be used for CMG momentum unloading. The effect of a constant aerodynamic 
torque on a gravitationally stabilized spacecraft was first studied by Garber [1]. 
Such torques produce constant attitude angles for which aerodynamic and gravi- 
tational torques are balanced. Garber has shown that small roll/yaw librational 
motions are affected by large pitch angles. The aerodynamic disturbance torques 
acting on the Space Station are expected to have constant values plus periodic com- 
ponents caused mostly by the effects of solar panel rotations and Earth’s diurnal 
bulge. As a result, attitude and CMG momentum oscillation about the torque equi- 
librium attitude will occur. A recent study [2] demonstrates the usefulness of the 
linear-quadratic-regulator synthesis technique and the concept of periodic distur- 
bance accommodation in minimizing attitude and/ or CMG momentum oscillations 
as needed for mission requirements. 

This paper is primarily concerned with attitude control and CMG periodic dis- 
turbance rejection for large-angle pitch maneuvers of the Space Station. New results 
and control concepts are used to extend the control scheme developed in [2]. Pitch- 
coupled roll/yaw equations of motion, first discussed in [l], and written in terms 
of Euler angles, are derived here in terms of quaternions. It is shown that these 
equations are well suited for use in designing a roll/yaw controller for large pitch 
motions of the Space Station. A simple concept of using quaternions for the con- 
trol of spacecraft large-angle maneuvers has been developed in [3, 4]. The concept 
is extended here to a more complicated case of controlling both the attitude and 
CMG momentum of the Space Station. Furthermore, this paper presents a new 
control concept of asymptotic momentum management of the CMGs, which tunes 
the open-loop unstable vehicle to a stable oscillatory motion during steady-state 
operations, while minimizing the control effort needed. 

Figure 1 is a functional block diagram representation of a quaternion feedback 
control system proposed for the Space Station. The attitude determination sys- 
tem utilizes rate gyros and star trackers to compute inertial quaternions and the 
absolute angular velocity of the Space Station. Relative quaternions with respect 
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to a local reference frame are then computed for control purposes. The proposed 
attitude/momentum controller utilizes relative quaternions, body rates (or relative 
quaternion rates), and CMG momenta to generate proper control torque commands 
to the CMGs. 


Mathematical Models 


In this section, equations of motion for the Space Station in a circular orbit are 
derived in terms of quaternions. For simplicity, the Space Station is assumed to 
be a single rigid body. Emphasis is on the use of quaternions in the equations of 
motion and in feedback control. Quaternions define the rigid body attitude as an 
Euler-axis rotation. The vector part of the quaternions indicates the direction of the 
Euler axis. The scalar part of the quaternions is related to the rotation angle about 
the Euler axis. Detailed discussion of the kinematics associated with quaternions 
and Euler angles can be found in many texts (see, e.g., [5, 6]). 

The relationships between quaternions and Euler angles, for the pitch-yaw-roll 
body-axis rotation sequence used in this paper, are 


7i 


72 


?3 


. ?4 . 



cos(#i/2) sin(0 2 /2) sin(0 3 /2) + sin(<?i/2) cos(0 2 /2) cos(0 3 /2) 
cos(0j/2) sin(0 2 /2) cos(0 3 /2) + sin($i/2) cos(0 2 /2) sin(0 3 /2) 
c°s(0i/2) cos(0 2 /2) sin(0 3 /2) — sin(0!/2) sin(0 2 /2) cos(0 3 /2) 


( 1 ) 


where (0j, $ 2 , 9 3 ) are the roll, pitch, and yaw Euler angles of the body axes with 
respect to the local vertical and local horizontal (LVLH) axes, which rotate with 
the orbital angular velocity; and (g 2 , q 2 , q 3 ) are the vector parts of the quaternions 
which indicate the direction of the Euler axis while q 4 is the scalar part of the 
quaternions and is related to the angle of rotation about the Euler axis. Inverse 
relations may also be written as: 


2(<7 i '74 — 7293) 


6\ = tan 1 

62 = tan -1 , 

[l-2ql-2ql 

63 = s\n~ l [2{ qi q 2 + q 3 q 4 )] . 


1 ~ 2<jj — 2ql J 
2(7274 — 7 i 7 a) 


( 2 ) 


The nonlinear equations of motion and attitude kinematics for the Space Station 
are as follows: : 
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Space Station Dynamics: 
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Attitude Kinematics (with respect to LVLH): 
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CMG Momentum: 
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( 4 ) 


( 5 ) 


and (wi, wi, u> 3 ) are the body-axis components of the absolute angular velocity of 
the Space Station; I, } (i = j) are the moments of inertia; I {j (i ^ j ) are the products 
of inertia; (h u h 2 , h 3 ) are the body-axis components of the CMG momentum; (u lt 
« 2 , « 3 ) are the body-axis components of the control torque; (w lf w 2 , w 3 ) are the 

body-axis components of the external disturbance torque; and n is the orbital rate 
of 0.0011 rad /sec. 

When body and control axes are aligned with the principal axes of the Space 
Station (/j = / U) J 2 = / 22) j 3 _ Eqs. ( 3 ) become 


Awi - (/ 2 - I 3 )u 2 u 3 + 6n 2 (/ 2 - I 3 )(qiq 4 + q 2 q 3 )(l - 2 q\ - 2 q\) = -u x + w, 

I 2 u 2 - ( I 3 - Ii)u)\U 3 + 6 ti 2 {I 3 - I\){q\ q 3 - q 2 q 4 )( 1 - 2q\ - 2 q\) - -u 2 4- w 2 (6) 
hu 3 - (J x - I 2 )uiu 2 + 12n 2 (/ 1 - I 2 ){qiq 3 - q 2 q 4 ){qiq 4 + q 2 q 3 ) = -u 3 + w 3 . 
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These equations can also be found in [ 5 ]. In this paper, however, we present a new 
set of equations determined by linearizing the above equations for the case of large 
pitch angles with small roll/yaw attitude changes. In this case, Eqs. (4), (5), and (6) 
can be linearized with respect to q\ and q 3 as follows: 

Space Station Dynamics: 

Ji<bi + n(/ 2 - I 3 )u> 3 + 6n 2 (/ 2 — I 3 )( 9 4 - 2q < q\)q l 

+ 6n 2 (J 2 - h){2q 2 q\ - q 2 )q 3 
/2W2 + 6n 2 (Ji - I 3 )(q 2 q 4 - 2 q A q%) = -u 2 + w 2 
^3^3 + n(Ii — /2V1 — 12n 2 (/x — I 2 ){q 2 q\)qi 

- 12 n 2 (/x - / 2 )(g49 2 )?3 = 

Attitude Kinematics (with respect to LVLH): 

. <7 4 ?2 

91 = + yw 3 + nq 3 

92 ~ + n ) 

■ 92 $4 

?3 = - y w i + y W 3 “ n?1 

94 = — “-(^2 + n ) 

CMG Momentum: 

hx — nh 3 = «x 

= u 2 (9) 

h 3 -f nhi = u 3 . 

For the case with large pitch angles and small roll/yaw attitude changes, the re- 
lationships between quaternions and Euler angles can be simplified by linearizing 
with respect to #1 and 6 3 . Equations (1) then reduce to 

• 02 

q 2 = sin y (10a) 

02 

q 4 = cos y (10 b) 
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Inverse relations for Eqs. (10c) are written as' 
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( 11 ) 


Equations (7) and (8) may be used to derive the Space Station equations of motion 
in terms of quaternions, written as follows: 
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The quaternion relations of Eqs. (10) may be used to transform Eqs. (12) to the 
following known form [1, 6] involving only Euler angles: 

I\0\ + n 2 (J 2 - / 3 )(1 + 3cos 2 0 2 )0i - «(/, - J 2 + I 3 )9 3 

+ 3 n 2 (J 2 - J 3 )(sin^ 2 cos^2)^3 = ~u x + W! (13a) 
I 2 O 2 + 3 n 2 (Ii - / 3 )sin0 2 cos0 2 = -u 2 + w 2 ( 13 ^ 

hh + n 2 (/ 2 - 7j)(l + 3sin 2 6> 2 )03 + n (I x - I 2 + I 3 )9 X 

+ 3 n 2 (/ 2 - I ! )(sin0 2 cos0 2 )0i = -«3 + w 3 . (13c) 

A final linearization with respect to small pitch motions leads to the following well- 
known equations of motion: 

Space Station Dynamics: 


I A + 4 n 2 (/ 2 - I 3 )9 l - n( I x - I 2 + I 3 )9 3 = - Ul + Wl 
I 2 9 2 + 3n (/1 — I 3 )9 2 = — u 2 + w 2 

+ n 2 (/ 2 - I x )9 3 + n(/j - I 2 + I 3 )9 X = -u 3 + w 3 

Attitude Kinematics (with respect to LVLH): 

9\ — n0 3 = 

9 2 — n = u 2 

^3 + n#! = o.’ 3 . 
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(14a) 

(14b) 

(14c) 


(15a) 

(15b) 

(15c) 


These linearized equations are used in [2] for the case of small roll, pitch, and 
yaw attitude changes. In this paper, emphasis is on the use of Eqs. (7) for the 
momentum/attitude control of the Space Station having small roll/yaw attitude 
changes but large-angle pitch motions. 

Inertia values for the Phase 1 Space Station, as well as assembly flight 3, are 
listed in Table 1. Table 2 includes expected aerodynamic disturbances which are 
modeled as a bias plus periodic terms in the body-fixed control axes: 

w(t) = Bias + A n sin(nt+<£i) + A 2n sm(2nt+£ 2 ) + A 3n sm(3nt+4 > 3 ) +A 4n sin(4nt+<j> 4 ) 

The disturbance torque acting in each axis is determined from data generated at 
NASA Johnson Space Center by a nonlinear simulation program. The program 
simulates translational and rotational motions of the Space Station in orbit about 
an oblate Earth. It includes rotating solar panels, time- varying surface areas, and 
time-varying center-of-pressure locations. A Jacchia-Lineberry atmospheric model 
.Itw compute density variations. The dominant aerodynamic torque frequen- 
cies at n and 2n are caused by Earth’s diurnal bulge and solar panel rotation 
effects, respectively. Actual magnitudes and phases of these disturbance torques 
are assumed unknown for control design. 

Control Issues 

Before presenting the pitch and roll/yaw controller designs, it is important to 
clarify some issues related to the effects of large pitch motions and inertia value 
uncertainties on the stability of the controlled Space Station. A characteristic of 
momentum/attitude control using gravity-gradient torque is that pitch, roll, and 
yaw responses will settle down to, or oscillate about, a constant torque equilibrium 
attitude (TEA). Primary factors involved in determining the constant pitch TEA 
are the magnitude of the bias in the disturbance torque and the numerical difference 
between roll and yaw moments of inertia. This can be seen by studying the steady- 
state form of Eq. (14b): 

6 = w 

2 3 n 2 (/, - h) 



w ere S, . S the pitch TEA angle and w is the bias of the pitch disturbance torque. 
Pitch gravity-gradient torque is largest when the pitch attitude is 45 degrees This 
is predicted in Eq. (13b) where, at the steady-state, 

2w 


sm 


(2<i)) = 


3"Vi - h) ' 

It may be necessary to consider this worst-case pitch TEA in control system design. 
After switching to CMG mode from some other modes (e-g„ reboost mode utilizing 
reaction jets), the Space Station must be able to achieve TEA in each axis without 
CMG momentum or commanded torque saturations. 

The large motions possible in the pitch axis emphasize the importance of using 
Eqs. (7), (12), or (13) in designing the control system. These equations show the 
dependence of roll/yaw dynamics on pitch attitude. In fact, a roll/yaw closed- 
loop system designed for small pitch angles may become unstable at large pitch 

angles; therefore, roll/yaw closed-loop stability must be checked at various pitch 
TEA values. 


Other factors affecting closed-loop stability are uncertainties and variations in 
moments and products of inertia. In particular, the magnitude of the pitch gravity, 
gradient torque depends on the difference between the roll and yaw moments of 
inertia. If these inertias are nearly equal (e.g., see Table 1, assembly flight 3 data), 
pitch gravity-gradient torque is very small and the pitch TEA (if it exists) is large. 
In addition, very small uncertainties in the moments of inertia can cause the system 
to become unstable. A discussion of the importance of checking closed-loop system 
obustness with respect to inertia uncertainties, by varying inertia values in an 
appropriate ‘direction’, is included in the appendix of this paper. 

Under normal operating conditions, the Space Station will have rotating solar 
y This causes time-varying (sinusoidal) roll and yaw moments of inertia, and 
consequently, a similarly time-varying gravity-gradient torque in the pitch axis as 
shown Fig. 2. If the bias value of /,( f)-/ 3 «) is such that the pitch gravity-gradient 
torque never changes sign, the system will remain stable. If the gravity-gradient 
torque does switch signs, the system may become unstable, depending upon the how 
long the sign of the torque is changed. For sufficiently short periods of this opposite 
torque, the system will be stable with large, bounded responses about the TEA. If 
gravity-gradient torque is to ultimately be used in momentum/attitude control, the 
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above factors will be useful in the design of the inertia configuration. Consideration 
should be given to defining a boundary near the point or condition corresponding 
to zero pitch-axis gravity-gradient torque, defined by I^t) - J 3 (<) = 0. This is 
illustrated in Fig. 2. The boundary may be thought of as the minimum gravity- 
gradient torque allowed in the pitch axis for which closed-loop stability, with respect 
to inertia uncertainties, is maintained. 

In the next sections, the pitch and roll/yaw controller designs are presented 
along with time simulations of the closed-loop system. In this paper, we expand 
on the previous study [2] by developing a periodic maneuvering controller which 
achieves asymptotic momentum management of the CMGs in all three axes. The 
Phase 1 inertia configuration listed in Table 1 is used. Effects of products of inertia 
are assumed negligible. The corresponding time-varying gravity-gradient torque 
in the pitch axis does not change sign; however, the time simulations presented 
here use the assumption that inertia values remain constant while periodic terms 
in the aerodynamic torque include solar panel rotation effects. Large pitch TEA 
responses are produced by introducing an appropriately large bias in the pitch-axis 
disturbance torque. 


Pitch Control 

In this section, the pitch-axis controller is developed for attitude and CMG 
momentum control. It is shown in [2] that disturbance rejection filters can be used 
to reject ei .her attitude or CMG momentum oscillations occurring at the frequencies 
present in the disturbance torques. Since asymptotic momentum management of 
the CMGs in all three axes is of primary interest in this paper, the disturbance 
rejection filters for the pitch axis have the following forms: 

d' 2 + (n) 2 a 2 = h 2 
02 + (2 Tl) 2 02 = 

72 + (3n) 2 7 2 = h 2 
+ (4n) 2 77 2 = h 2 

where initial conditions for the filter states can be arbitrarily selected (usually zero 
initial conditions). Use of filters at frequencies n, 2n, 3n, and 4n is indicated by 
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aerodynamic torque data generated by a nonlinear simulation program written for 
the Space Station. The pitch-axis control logic is given by a single control input 
involving twelve states: 

u 2 = K 2 x 2 (16) 

where 

_ _ A 

K 2 = a 1 x 12 gain matrix 

x 2 = [ ?2 92 h 2 f h 2 a 2 a 2 fi 2 /3 2 j 2 j 2 q 2 rj 2 ] T . 

The control task is to find proper gains for this twelve-state feedback controller. 

In order to use linear control design methodologies, Eq. (7b) must be linearized 
for small pitch motion. This results in Eq. (14b), which is used as the basisTor pitch 
control analysis and design. Various techniques may be used in selecting the twelve 
gains of Eq. (16). These include linear-quadratic-regulator (LQR) synthesis [7] 
and direct assignment of closed-loop eigenvalues using a pole-placement technique. 
Several iterations of any method may be required to achieve satisfactory closed-loop 
performance and robustness. Note that gains resulting from Eqs. (14b) are for the 
use of Q 2 in state feedback. In order to accommodate q 2 for use in feedback, the 
gains corresponding to states q 2 and q 2 are doubled since the approximation used 
for q 2 is & 2 j 2. New gains do not need to be computed for the pitch controller in 
the case of large pitch motion. It is mostly roll/yaw destabilization at large pitch 
angles which forces pitch-axis instability. 

The open-loop pitch axis of the Phase 1 inertia configuration is unstable, with 
poles at s = ±1.5n, 0, 0, and filter poles at s = ±jn, ±j2n, ±j3n, ±j4n. One pole 
at s = 0 comes from the integral feedback of h 2 . After iterative use of an LQR 
synthesis code, available in CTRL-C software, a set of closed-loop eigenvalues have 
been selected and are listed in Table 3. The corresponding gain set is given in 
Table 4. Closed-loop pitch responses of Eq. (6), for a pitch-axis maneuver of —30 
degrees (caused by a large pitch-axis torque bias), can be seen in Figs. 4, 5, and 6. 
Comments on the responses are reserved until after the presentation of the roll/yaw 
controller design. 


Roll/ Yaw Control 

The roll/yaw controller has a structure similar to that of the pitch controller. By 
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examining the open-loop transfer function matrix from control inputs to roll/yaw 
attitudes and CMG momentum, it is shown in [2] that a periodic disturbance at the 
orbital rate.can be rejected in the yaw attitude but not in the roll attitude. The 
analysis is accomplished using Eqs. (9) and (14), which assume small motions in all 
axes. In an effort to determine if it is possible to have periodic-disturbance rejection 
in both roll and yaw attitudes for this different case involving large pitch motions, 
a similar analysis is considered here using the pitch-coupled roll/yaw dynamics 
described by Eqs. (13). 

By combining Eqs. (9) and (13), the transfer function matrix from (u x ,tX 3 ) to 
(0i,& 3 ) can be written as 


' *»(«) ' 

1 


~ A 




G 13(5) 
^33(5) 


u x (s) 

u 3 (s) 


where 


Gn(s) = -[J 3 s 2 + (1 + 3s 2 )n 2 (J 2 - J x )][ 5 2 + n 2 ] 

Giz(s) = — [n(/j — J 2 + h) s — 3n 2 (/j — / 3 )s 2 c 2 ](s 2 + n 2 ] 

G 3 i(s) = [n(/ x — I 2 + l 3)3 + 3n 2 (/ 2 — Ji)s 2 c 2 ][s 2 -f- n 2 ] 

<^ 33 ( 5 ) = -[hs 2 + (1 + 3c 2 )n 2 (I 2 - / 3 )][s 2 + n 2 ] 

A = A J 3 O s 2 + n 2 ){s 4 + n 2 [l + 3* x + ^*3 + 3 (k 3 - Ms 2 ]* 2 + 

[3n 3 (& 3 — &i)s 2 c 2 ]s + 4n 4 fc x & 3 } 

and s 2 = sin# 2 ,c 2 = cos# 2 ,fc x = (J 2 — i 3 )// x , and k 3 = (J 2 — / x )// 3 . Transmission 
zeros at ±jn appear in the transfer function matrix. It would seem that periodic 
disturbances of frequency n cannot be rejected in either the roll or yaw attitudes; 
however, for CMG momentum and control torque relations defined by 

A x = u x 
A 3 = «3 


and with appropriate alterations of Eqs. (13), the transfer function matrix from 
( u i> u 3 ) to (0 X ,# 3 ) can be written as 


' 9M ' 

1 


~ A 


G u (s) 

G 3 l (s) 


Gi 3 (s) 
G33 (s) 


u x (s) 

“ 3 ( 5 ) 
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where 


G u = 0){-/ 3 s 3 + n 2 [3(/i - I 2 )s\ - I 3 ]s + 3n 3 (/ 2 - / 3 ) 52 c 2 } 

Gn = (ns){(/ 2 - h)s 2 + [3n(I 2 - J 3 )s 2 c 2 ]s + (1 + 3 s 2 )n 2 (J 2 - /j)} 

G 3l = — (ns){(/ 2 - / 3 )s 2 + (1 + 3c 2 )n 2 (/ 2 - / 3 )} 

G 33 = -{hs* + n 2 [3(/ 2 - I 3 )c 2 + h]s + 3 n 2 (J 2 - A> 2 c 2 > 

A = J,/ 3 (s 2 ){s 4 + n 2 [ 1 + 3^i + *^3 + 3(Jfc 3 - *x)s 2 ]s 2 + 

[3n 3 (k 3 - ^> 2 c 2 ]s + 4n 4 M: 3 }. 

Transmission zeros are not apparent in these expressions. A numerical analysis 
reveals, however, that there are transmission zeros at ±jn for the transfer functions 
from «i and u 3 to roll attitude, while yaw attitude has no troublesome zeros. These 
results show that, even for the case of pitch-coupled roll/yaw dynamics, there is 
an inability to reject roll attitude oscillations occurring at the orbital rate. Hence, 
in this paper, periodic-disturbance rejection for CMG momentum in both the roll 
and yaw axes is considered. That is, a periodic control of the Space Station for 

asymptotic momentum management of the CMGs in all three axes is of primary 
interest here. 

Periodic-disturbance rejection filters for the roll/yaw axes can be represented 
as: 


&i + (n) 2 oti = hi 
+ (2 n) 2 0i = hi 
7i + (3n) 2 7j = hi 
*7i + (4 n) 2 rji = hi 
&3 + (n) 2 a 3 = h 3 
k + (2 n) 2 0 3 = h 3 
73 + (3n) 2 7 3 = h 3 
V3 + (4n) 2 77 3 = h 3 . 


The roll/yaw control logic involving two control inputs and twenty-four states 
is expressed as 


Ui 


= I< 


*1 

X 3 


( 17 ) 
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where 


0i 7i 7i 9i i?i j (roll states) 

03 73 73 *13 V 3 J (yaw states) . 

The pitch-coupled roll/yaw equations described by Eqs. (7) are used as the basis 
for the roll/yaw controller design. It is evident in the pitch-coupled equations that 
roll/ yaw dynamics are dependent upon pitch attitude. Consideration must be given 
to the pitch-axis TEA when designing the roll/yaw controller gains. Equations (7) 
(or Eqs. (12) or (13)) are especially useful for this purpose when q 2 and q 4 in the 
equations are assigned their respective values corresponding to the expected pitch 
TEA (see Eqs. (10a) and (10b)). As a result, Eqs. (7a) and (7c) become linear 
and any linear control design methodologies may be used to design the roll/yaw 
controller gains. 

For spacecraft influenced by gravity-gradient torques, it is interesting to examine 
the changes in roll/yaw open-loop eigenvalues that occur as the pitch bias changes. 
This was first studied by Garber [l]. A root locus of open-loop eigenvalues versus 
pitch angle, for the Phase 1 inertia configuration, is shown in Fig. 3. It can be seen 
that the open-loop roll/yaw dynamics are not very sensitive to pitch attitude. The 
Space Station is unstable with poles at s = ±1.05n±j0.7n, 0, 0, ±jn, and filter poles 
a ! s = ±j2n, ±j2n, ±j3n, ±j3n, ±j4n, ±j4n (for 8 2 = 0° where g 2 =0 and 

$ 4 = 1 ). The double pole at s = 0 occur because of the integral feedback of hi and 

- -Vi, * -• t X i - - - - 

After iterative use of an LQR synthesis code, closed-loop eigenvalues have been 
selected and are listed in Table. 3. A gain set for Eqs. (17), corresponding to a 
pitch TEA of 0°, is listed in Table. 4. For these gains, the closed-loop roll/yaw 
axes are stable for pitch angles ranging from -21° to +23°. Since a simulation of 
the large-angle pitch maneuver needed to reach a pitch TEA of -30° is desired, 
a different gain set is used for the simulations presented in this paper. For these 
gains, the closed-loop roll/yaw axes are stable for pitch angles ranging from -4S° 
to +3°. Closed-loop roll/yaw responses of Eqs. (6), for a pitch-axis maneuver of 
—30 degrees, can be seen in Figs. 4, 5, and 6. The overall closed-loop system has a 
10 dB gain margin and a phase margin of 60° in each control loop. 
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K = a 2 x 24 gain matrix 

Xi = [ qi u>! hi f hi at t di 

x 3 = [ 93 h 3 f h 3 q 3 d 3 /? 3 


Discussion of Simulation Results 


Closed-loop responses for a simulation of the nonlinear dynamics described by 
Eqs. (6) are shown in Figs. 4, 5, and, 6. Quantities plotted include quaternions, 
CMG momenta, and control torques. Initial conditions corresponding to 0 t (O) = 
0 2 (O) = *(0) = 1 degree and 0,(0) = 0 2 (O) = 0 3 (O) = 0.001 deg/sec are assumed. 
Allowable limits on CMG momentum and commanded torque are assumed to be 
30,000 ft-lb-sec and 150 ft-lb, respectively. 

In the roll axis, quaternion oscillates (±0.7°) about a roll TEA of -0.003 
(0, = -0.5°). Roll CMG momentum h, is the input to the roll-axis disturbance 
rejection filter and settles down to zero after reaching a maximum value near 8000 
ft-lb-sec. Control torque u, is zero at the steady-state and has a peak value near 
23 ft-lb. In the pitch axis, quaternion q 2 oscillates (±4.3°) about a pitch TEA of 
-0.257 (0 2 = -30°). The large-angle maneuver causes the pitch CMG momentum 
h 2 (disturbance filter input) to become quite large at nearly 18,000 ft-lb-sec, before 
settling to zero. Control torque u 2 is zero at the steady-state with a maximum value 
near 30 ft-lb. In the yaw axis, quaternion q 3 oscillates (±1.2°) about a yaw TEA 
of —0.013 (0 3 = -1.5°). Yaw CMG momentum h 3 (disturbance filter Input) settles 
down to zero after reaching a maximum value close to 600 ft-lb-sec. Control torque 
u 3 is zero at the steady-state, and reaches a maximum value near 12 ft-lb. 

The simulations show that the proposed control scheme tunes the open-loop 
unstable Space Station to a stable, oscillatory motion which minimizes control effort 
aur ng steady-state operations. For the assumed disturbance torque models (with 
unknown magnitudes and phases), the stabilized Space Station needs no control 
torque at steady-state conditions. Analysis shows, however, that small-amplitude 
periodic components of frequencies 5n and 6n are present in u, and u 3 at the steady- 
state. These small residual components are caused by the coupling between the pitch 
and roll/yaw axes, and become particularly noticeable for large pitch biases. 

There is an interesting feature of the quaternion feedback scheme which is not 
apparent from the simulation responses. If rejection of pitch and yaw attitude 
oscillations is desired, it seems natural to use q 2 and q 3 as inputs to the respective 
disturbance rejection filters. Even though q 2 and q 3 will become constant at the 
steady-state, all of the Euler angles will oscillate. A study of Eqs. (2) (or Eqs. (11)) 
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reveals why. By assigning constant values to q 2f ? 3 , and q 4 in Eqs. (2) (or Eqs. (11)), 
it can be seen that 9 U 0 2 , and 0 3 are ail functions of q x , which oscillates. The same 
may be said for q u q 2 , q 3 , and if Euler angle feedback is used. For oscillations of 6 X 
and constant values for 0 2 and 0 3 , Eqs. (1) show that all quaternions are functions of 
0i and will therefore oscillate. In either case, however, these oscillations are small. 
The important point is that the elimination of pitch or yaw oscillations, if needed, 
may be accomplished by using Euler angles 0 2 and 0 3 as disturbance rejection filter 
inputs (with appropriate gain changes). 

A check of closed-loop robustness with respect to inertia uncertainties empha- 
sizes an important issue associated with the Phase 1 inertia configuration. By se- 
lecting specific ‘directions’ in which to vary the three moments of inertia (J X) / 2 , J 3 ), 
the closed-loop system can be shown to be unstable for as little as -7% uncertainty 
in / 3 with +8% uncertainty in I x . For these inertia variations, the pitch-axis gravity- 
gradient torque disappears (I x - 7 3 = 0) and closed-loop pitch dynamics become 
unstable. The limitations shown in this example (and several others involving even 
smaller inertia uncertainties), are not related to the selection of control logic but 
are physical limitations inherent to the current inertia configuration of the Phase 
1 Space Station. A description of the inertia variation ‘directions’ used above is 
presented in the appendix of this paper. 

Conclusions 

This paper has presented a new control concept for the planned Space Station 
involving asymptotic momentum management of control moment gyros. It has 
been shown that the proposed controller tunes the Space Station, which has a 
gravitationally ‘unstable’ inertia configuration, to a stable, oscillatory motion which 
minimizes the control effort needed at the steady-state. By utilizing the concepts of 
quaternion feedback control and periodic-disturbance rejection filters, the proposed 
controller provides robust control of the Space Station for large-angle pitch motions. 
The pitch-coupled roll/yaw equations of motion derived in this paper have been 
shown to be particularly useful in roll/yaw controller design and stability analysis. 


Appendix: Inertia Variations for Checking Closed-Loop Stability 

It is a common practice in control design to satisfy time and frequency-domain 
requirements first, then check for closed-loop robustness. For spacecraft, inertial 
properties may be very sensitive parameters in the closed-loop system. It seems 
reasonable to increase and decrease all inertias by the same percentage, thus check- 
ing controller effectiveness for a proportionally heavier and lighter spacecraft. This 
procedure may not indicate the true inertia sensitivity. It is important to consider 
the magnitude and direction of the variation for each inertia value. Since the inertia 
matrix may be transformed to three principal moments of inertia by aligning the 
body and principal axes, suggestions for variations in the moments of inertia for 
the roll (Ij), pitch ( I 2 ), and yaw (/ 3 ) axes are presented here. 

Three important relationships may be derived from the definitions for the mo- 
ments of inertia. These relationships axe as follows: 

A Ij > Iz) I\ + I3 > I2 , I2 + I3 > Jj . 

Together, these relations define the physically possible inertia configurations. A 
control designer may unknowingly use inertia venations which result in inertia val- 
ues that violate these constraints. Stability of the closed-loop system will be tested 
for a physically impossible inertia configuration. The important point is not the 
fictitious inertias but whether or not the control designer can redirect this extrane- 
ous stability margin to encompass more of the region of physically possible inertia 
values. 

When gravity-gradient torque is used in the control of a spacecraft, additional 
inertia constraints are introduced. The control scheme presented in this paper is a 
good example. Equations (14) show that the additional constraints are I 2 ^ I 3 and 

h ± h since I 2 = I 3 and 7 X = I 3 correspond to zero gravity-gradient torque in the 
roll and pitch axes, respectively. 

A useful aid for visualizing the relationship between inertia constraints and iner- 
tia variations is now presented. Figure A.l shows a three-dimensional, cubic figure 
defined in three ‘inertial’ directions. The inertia constraint relations may now be 
visualized as planes in this ‘inertial’ space. The planes I x + / 2 = / 3i /, -f-/ 3 = / 2> and 
Ii + h = h are labeled in Fig. A.l, and represent the physical boundaries of inertia 
values. The area inside the three intersecting boundaries represents all physically 
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possible inertia configurations. A representation of the cut-away portion of the ‘in- 
.^rtMLeube’ is shown in Fig. A.2. Planes defining the physical and gravity-gradient 
boundaries are labeled. Spacecraft inertia configurations which are open-loop stable 
for gravity-gradient control schemes may be identified by requiring positive coeffi- 
cients for the open- loop pitch, and roll/yaw characteristic equations (see Fig. A. 2). 
The resulting stable configurations are represented by the region in Figs. A.l and 
A.2 where I 2 > I\ > / 3 . Figures A.l and A. 2 provide a three-dimensional represen- 
tation of the information presented in the fc 3 versus k A inertia ratio plots found in 
[ 5, 6}. It may be convenient to normalize the moments of inertia being studied by 
\/^i + -^2 + -I? i n order to locate the position of the nominal configuration within 
a ‘unit inertia cube’. The relative positioning of the nominal inertia configuration 
from the constraint boundaries can then be easily determined. 

Since the shortest distance from a point to a plane is in a direction perpendicular 
to that plane, it seems logical to check inertia variations in directions perpendicular 
to the inertia boundaries. In this way, the minimum variation necessary to reach 
a physical boundary can be found while checking the closed-loop stability of the 
system in question. For gravity-gradient control systems, inertia variations perpen- 
dicular to the planes I\ = I 2 and I 2 = J 3 are needed. It may be seen in Fig. A.2 
that these planes intersect inside the region of physically possible inertia values, 
and partition the region into several sections. It should be a control designers’ goal 
to include the area within the physical boundaries inside a ‘control surface’ which 
contains all of the inertia values for which the closed-loop system is stable. 
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Table 1. Space station inertia 
configurations 


Inertia 

(slug-ft 2 ) 

Assembly 

Flight#3 

Phase 1 

/ ii 

23.22E6 

50.28E6 

27 

1.30E6 

10.S0E6 

l 33 

23.23E6 

58.57E6 

In 

— 0.023E6 

-0.39E6 

hz 

0.477E6 

0.16E6 


— 0.011E6 

0.16E6 


Table 2. Phase 1 aerodynamic torque models 
(in units of ft-lbs) 

w i l-t-sin(nt)-f0.5sin(2nt) 

+O.3sin(3nt)+0.5sin(4nt) 
w 2 13* + 1.2sin(nt)+3.5sin(2nt) 

-f-0.3sin(3nt)+0.5sin(4nt) 
w 3 l+sin(nt)+0.5sin(2nt) 

-f 0.3sin(3nt)+0.5sin(4nt) 

* nominal pitch bia_s torque is 4 but 13 is used 
to produce a large pitch TEA 


Table 3. Phase 1 closed-loop eigenvalues 
(in units of orbital rate - 0.0011 rad/sec) 



Momentum/Attitude 

Disturbance Filters 

Pitch 

-1.0, -1.5 

— 1.5±jl.5 

— 0.3±jl.0 
— 0.3±j3.0 

— 0.3±j2.0 
— 0.3±j4.0 

Roll/ Yaw 

-0.23,-0.71 
— 1.04±j0.70 

— 0.53±jl.54 
— 1.06±j0.71 

— 0.14±j0.99 

— 0.19±j2.01 

— 0.32±j3.02 
— 0.53±j3.97 

— 1.13±j0.75 

— 0.47±j2.20 

— 0.6S±j3.21 
— 0.25±j4.00 




Table 4. Phase 1 controller gains 


Pitch 

Units 

Roll/ Yaw* 

[K,f 


[KF 

4.2425E+2 

(ft-lb/rad) 

3.8526E+3 

3.7381E+2 

2.5412E+5 

(ft-lb-sec/rad) 

1.2003E+6 

1.0126E+5 

1.4840E— 2 

(ft-lb/ft-lb-sec) 

1.4360E— 2 

1.9364E— 3 

4.0150E— 6 

(ft-lb/ft-lb-sec 2 ) 

— 1.6361E— 6 

2.7852E— 7 

— 1.9064E— 9 

(ft-lb-rad 2 /ft-lb-sec 3 ) 

3.6578E— 10 

— 1.8526E— 10 

2.1970E— 6 

(ft-lb-rad 2 /ft-lb-sec 2 ) 

7.6282E— 7 

1.1857E— 7 

— 4.6097E— 9 

(ft-lb-rad 2 /ft-lb-sec 3 ) 

— 3.2712E— 9 

— 5.7517E— 10 

— 5.2383E— 7 

(ft-lb-rad 2 /ft-lb-sec 2 ) 

— 3.3865E— 7 

— 1.6409E— 7 

— 5.3793E— 9 

(ft-lb-rad 2 /ft-lb-sec 3 ) 

— 1.0702E— 8 

— 1.1317E— 9 

— 1.9423E— 6 

( ft-lb-rad 2 / ft-lb-sec 2 ) 

— 3.4827E— 6 

— 5.3664E— 7 

— 7.3458E— 9 

(ft-lb-rad 2 /ft-lb-sec 3 ) 

— 1.5903E— 8 

— 1.5491E— 9 

— 2.6056E— 6 

(ft-lb-rad 2 / ft-lb-sec 2 ) 

— 3.1256E— 6 

— 4.7197E— 7 



9.4016E+2 

2.4994E+2 



-1.2743E+5 

1.1386E+5 



— 2.4992E-3 

— 3.5209E— 3 



-7.3398E-7 

— 1.0348E— 6 



4.8557E— 9 

— 5.5935E— 10 



3.7017E— 7 

— 4.2651E— 6 



2.0608E— 9 

— 6.8224E— 10 



1.8854E— 6 

— 2.4769E— 6 



4.0142E— 10 

9.4962E— 10 



1.5548E— 6 

— 2.7820E— 6 



8.3363E— 10 

8.3453E— 10 



1.3125E-6 

— 2.5757E— 6 


* Designed for a pitch TEA of 0° 


Units 


(ft-lb/rad) 

(ft-lb-sec/rad) 

(ft-lb/ft-lb-sec) 

(ft-lb/ft-lb-sec 2 ) 

(ft-lb-rad 2 / ft-lb-sec 3 ) 

(ft-lb-rad 2 /ft-lb-sec 2 ) 

(ft-lb-rad 2 /ft-lb-sec 3 ) 

(ft-lb-rad 2 /ft-lb -sec 2 ) 

(ft-lb-rad 2 /ft-lb-sec 3 ) 

(ft-lb-rad 2 /ft- lb-sec 2 ) 

( ft-lb-rad 2 /ft-lb-sec 3 ) 

(ft-lb-rad 2 /ft-lb-sec 2 ) 

(ft-lb/rad) 

(ft-lb-sec/rad) 

(ft-lb/ft-lb-sec) 

(ft-lb/ft-lb-sec 2 ) 

(ft-lb-rad 2 / ft-lb-sec 3 ) 

( ft-lb-rad 2 /ft-lb-sec 2 ) 

(ft-lb-rad 2 /ft-lb-sec 3 ) 

(ft-lb-rad 2 /ft-lb-sec 2 ) 

(ft -lb rad 2 /ft-lb-sec 3 ) 

(ft-lb-rad 2 /ft-lb-sec 2 ) 

(ft-lb-rad 2 /ft-lb -sec 3 ) 

(ft-lb-rad 2 /ft-lb-sec 2 ) 



B/I indicates body frame with respect to inertial frame 
B/L indicates body frame with respect to LVLH frame 
I/L indicates inertial frame with respect to LVLH frame 


Fig. 1. A quaternion feedback control system 

for the Space Station. 































Open-loop characteristic equations for a spacecraft with 
gravity-gradient control scheme: 


Pitch axis: 


Roll/yaw axes: 


h{s 2 + 3 n 2 k 2 } = 0 


Iih{s* 4- n 2 [l + 3 ki + kikz]s 2 + 4n A k\k%} = 0 
where k t = (J 2 - I 3 )/I u k 2 = (h - / 3 )// 2 , and k 3 = (I 2 - 


+ ^2 — -^3 


h>h> li 


J3 > -?2 - M . 


h>I 3 > I , 


I 2 > h> I t 


^ \A > h > h h > ^ 


/1 = I 2 

/ 3 = o 


Fig. A. 2. Regions within the physical boundaries. 





